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We present a simple quantum many-body system - a two-dimensional lattice of qubits with a 
Hamiltonian composed of nearest-neighbor two-body interactions - such that the ground state is a 
universal resource for quantum computation using single-qubit measurements. This ground state 
approximates a cluster state that is encoded into a larger number of physical qubits. The Hamiltonian 
we use is motivated by the projected entangled pair states, which provide a transparent mechanism 
to produce such approximate encoded cluster states on square or other lattice structures (as well as a 
variety of other quantum states) as the ground state. We show that the error in this approximation 
takes the form of independent errors on bonds occurring with a fixed probability. The energy gap 
of such a system, which in part determines its usefulness for quantum computation, is shown to be 
independent of the size of the lattice. In addition, we show that the scaling of this energy gap in 
terms of the coupling constants of the Hamiltonian is directly determined by the lattice geometry. 
As a result, the approximate encoded cluster state obtained on a hexagonal lattice (a resource that 
is also universal for quantum computation) can be shown to have a larger energy gap than one on 
a square lattice with an equivalent Hamiltonian. 

PACS numbers: 03.67.Lx 



The state of a quantum many-body system can serve as 
a universal resource for quantum computing, where com- 
putation proceeds through single-qubit measurements 
alone This observation raises the intriguing pos- 

sibility that there exist physical systems which are "nat- 
urally" quantum computers. More precisely, one may ask 
whether there exist quantum many-body systems which, 
when cooled sufficiently close to the ground state, can be 
used for quantum computation by simply making indi- 
vidual measurements on the constituent particles. 

It is straightforward to write down a Hamiltonian for 
a many-body spin system for which this is the case - one 
needs only to take the (negative) sum of the stabilizer 
operators corresponding to, say, a square lattice cluster 
state 0. This Hamiltonian is somewhat unsatisfactory, 
however, because it involves five-body interactions, and 
fundamental interactions are strictly between_two bod- 
ies. Using the results of Haselgrove et al. Nielsen 
proved the following negative result: a cluster state suit- 
able for quantum computation cannot arise as the exact 
ground state of any Hamiltonian involving only local two- 
body interactions However, Oliveira and Terhal 9, 
building on the results of Kempe, Kitaev and Regev 3, 
demonstrated that cluster states (and other such states 
that are universal) can be approximated by the ground 
state of a local two-body Hamiltonian. The key idea of 
their result is to make use of "mediating" ancilla qubits to 
create an effective many-body coupling out of two-body 
interactions. 

Here, we provide a simple and explicit scheme for 
cluster-state quantum computation using the ground 
state of a Hamiltonian consisting of only two-body inter- 



actions. Our approach, which provides an alternative to 
the use of mediating systems as in 0, 0] , is motivated by 
the projected entangled pair states (PEPS) as proposed 
by Verstraete and Cirac '^l. Using PEPS, it was demon- 
strated that the cluster state can be obtained through 
a projection of a number of virtual qubits prepared in 
maximally- entangled states down to a lesser number of 
physical qubits. We turn this construction around, and 
consider how, through the use of nearest-neighbor two- 
body Hamiltonians, we may effect a similar projection - 
but in this instance of a number of physical qubits down 
to a lesser number of logical qubits. The system we in- 
vestigate is gapped, as the energy difference between the 
ground and first-excited states is independent of the size 
of the lattice. The low-energy theory is described by an 
effective five-body interaction Hamiltonian that is pre- 
cisely equal to the (negative) sum of the stabilizer opera- 
tors acting on the logical qubits. That is, a cluster state 
encoded as a logical state is approximated by the ground 
state of a larger number of physical qubits. This encod- 
ing provides a very transparent mechanism for achieving 
the effective many-body coupling out of two-body inter- 
actions, which can be of Ising and Heisenberg form, and is 
naturally extendible to other PEPS constructions. Cru- 
cially, despite using logical qubits consisting of several 
physical qubits, measurement-based quantum computa- 
tion can proceed on the encoded cluster state using only 
single-qubii measurements. 

We also investigate the usefulness of this ground state 
for measurement-based quantum computation. Specifi- 
cally, we analyze the errors associated with the fact that 
the ground state only approximates the encoded cluster 
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state, and demonstrate that the (fixed) energy gap is di- 
rectly determined by the lattice geometry. 

The PEPS on a square lattice Q consists of a 2- 
dimensional lattice of virtual qubits following the Archi- 
median tiling 4.8.8 (also known as a "CaVO" lattice) 0, 
as shown in Fig. ^ Quadruples of virtual qubits on the 
vertices of the 4- vertex tiles (such as those circled) form 
sites. Pairs of virtual qubits at neighboring sites, con- 
nected by dashed-line bonds in Fig. ^ are prepared in 
the two-qubit cluster state 
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The four virtual qubits at each site are projected down 
to a single physical qubit using the projector 



P = |0)(0000| + |i)(iiii| 



(2) 



The resulting state of the physical qubits on the lattice 
is the cluster state 

In our approach, we choose a nearest-neighbor two- 
body Ising-type Hamiltonian such as to effectively imple- 
ment the above projection on four physical, rather than 
virtual, qubits. Consider a lattice of physical qubits with 
the same lattice structure as the virtual qubits above, 
as in Fig. We label sites of four physical qubits by a 
Greek character, e.g., /i, and define S to be the set of all 
sites. Qubits at each site interact with nearest neighbors 
via a site Hamiltonian Hs which is of Ising form 
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where the first sum is over all sites, i ^ i' denotes pairs of 
neighboring qubits at site /i g S" connected by solid-line 
bonds in Fig.^ and erf is the Pauli z-operator for the 
physical qubit i at site /i. 

The bonds between sites in Fig. ^ indicate a different 
two-body interaction, given by a two-body Hamiltonian 
of the form 
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where is a physical Pauli x-operator, and (/i, i) ^ 
(h'jj) denotes pairs of neighboring qubits connected by 
dashed-line bonds in Fig. ^ 

The total Hamiltonian for the lattice is given by the 
sum of these terms, H = gHs + XV, where g and A 
have units of energy. To investigate the spectrum of this 
Hamiltonian, we will use perturbation theory in A/ g, with 
gHs the unperturbed Hamiltonian. 

We first investigate the spectrum of the unperturbed 
Hamiltonian gHs- At a single site, there are three energy 
levels. The ground-state is degenerate, two-dimensional, 
and spanned by the states 
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FIG. 1: Schematic of the spin lattice, with couplings. Solid 
lines represent an Ising-type interaction, as in the Hamilto- 
nian of Eq. The dashed lines represent a coupling as in 
the Hamiltonian of Eq. Groups of four physical qubits 
bound by the Ising interaction, as within the circle, form a 
site - a logical qubit. 



The ground state space of the unperturbed Hamiltonian, 
then, can be viewed as a logical qubit. The energy of 
this ground state space is —4(7. The first excited state is 
12-fold degenerate, and has a energy of 0. The second 
excited state is 2-fold degenerate and has a energy of Ag. 

Consider a lattice as in Fig. consisting of Ns sites 
with periodic boundary conditions. The spectrum of the 
unperturbed Hamiltonian gHs is straightforward. The 
ground-state space has energy E^^ = —AgNs, is 2^^- 
dimcnsional, and is spanned by all logical states of Ns 
qubits. We denote this space Hl- A convenient basis for 
this ground state space is as follows: let |C) denote the 
logical cluster state, and let \C{a,(3, . . .}) denote the logi- 
cal cluster state with logical Z-errors on the sites a, /3, . . .. 
(The logical Z operator is Z = |OL)(Oi| - |li)(li|.) The 
set of states {|C), |C{a}), |C{a, /?}),.. .}, running over 
logical Z-errors at all possible sites, forms a basis for this 
ground state space. 

The first-excited space is {12Ns ■ 2^^^^)-dimensional, 
and has energy e[^'' = —Ag{Ns — 1). Thus, for the un- 
perturbed Hamiltonian gHs, the gap from the ground to 
first-excited space is Ag. The second-excited space has 
energy £'2°'' — —4:g{Ns — 2). These energies will serve as 
the zeroth-order energies in perturbation theory for the 
total Hamiltonian. 

We now turn to perturbation theory. The term XV in 
the Hamiltonian, representing the coupling on the bonds, 
breaks the degeneracy of the ground state space; however, 
as we will see, it only does so at a fourth-order pertur- 
bation in X/g. Let H^ be the projection onto the ground 
state space of the unperturbed Hamiltonian gHs, i.e.. 
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onto the "logical" space. It is straightforward to show 
that IlLyilL = 0, and thus there is no first-order correc- 
tion to the energies. It is also straightforward to show 
that 

UlV'^Ul = 2 X (number of bonds) = 4NsUl . (6) 

(This constant arises from the fact that each Pauli term 
in V squares to the identity.) Thus, there is a constant 
second-order correction to the ground-state energy - a 
shift - given by 
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There is no third-order correction to the ground state en- 
ergy, following from the fact that HlV^Hl — 0. Finally, 
at fourth-order, the degeneracy is broken. To calculate 
the correction to the energies, we will make use of two 
convenient properties of our encoding. First, we note 
that the product of all four operators on a single site 
(one on each physical qubit) results in a logical X oper- 
ator 
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Second, we note that a single operator acting on any 
of the four qubits at a site is equivalent to a logical Z at 
that site. With these facts, we find that 

^lV^^l - {J^lV'^^l? = 4! ^ i^^ , (9) 



where 



(10) 



is a stabilizer of the logical cluster state. In this expres- 
sion, X^ is a logical X acting at the site /i, and Z^, is 
a logical Z acting on a site that is connected to /i by a 
bond. 

The operator 11^1^ ''lli is already diagonal in our cho- 
sen basis (the states of the form \C{a, (3, . . .})). The 
fourth-order energy corrections for these states are de- 
termined by their eigenvalues, which are straightforward 
to calculate using the properties of stabilizers. First, we 
note that the cluster state |C) is an eigenstate of all sta- 
bilizers in the sum Eq. ^ with eigenvalue -t-1, and thus 
the fourth-order correction for the energy associated with 
this state is 
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Next, consider a state |C{a}) = Za\C), a cluster state 
with a single Z-error at the site a. This state is also an 
eigenstate of all stabilizers in the sum Eq. @ with eigen- 
value -|-1 except the stabilizer Kq, which has eigenvalue 



— 1. (This result is due to the fact that X^ and Za an- 
ticommute.) Thus, to fourth order, the spectrum of the 
Hamiltonian H = gHg + XV is as follows: the component 
of the ground state in Hl is the cluster state |C), with 
energy 
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The nth excited space is (^^)-dimensional. The compo- 
nent of this space in Ti^ is spanned by states obtained 
from \C) with n logical Z errors. These states have en- 
ergy En = Eo + nA, up to n = Ns, where 



A = Ei-Eo 
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(Above n = Ns, there is a large gap of order 6g to the first 
manifold of illogical states, which are those that include 
a single physical error.) We note that the gap A is 
independent of the size of the lattice. Intuitively, then, 
one may associate logical Z errors on any site with a fixed 
energy A = 3X'^/{16g^) each. 

Although the ground-state energy degeneracy of the 
unperturbed Hamiltonian gHs is not broken until fourth 
order, the lowest energy eigenstates are corrected even 
at first order as they will include terms from the higher- 
energy eigenspaces of gHs- A simple counting of the 
states involved, and the magnitude of the components, 
will allow us to calculate the overlap of the first-order- 
corrected ground state with the exact cluster state. It 
will also allow for a simple determination of an error 
probability per bond associated with this correction. 

At first order, a state in the ground-state space (we will 
use |C) as an example, but this calculation holds for any 
state) will be corrected to include components from states 
in the first excited space of gHs whenever (?/;|y|C) is 
non-zero. We note that every term of the form or 

(g)a^ from V acting across a bond will connect |C) to a 
unique excited state, all of which form an orthogonal set; 
thus, summing over all the terms in V, there will be ANs 
orthogonal states, denoted |fc), for which (fc|l/|C) ^ 0. 
The first order corrected ground state \Eq) is 
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We now consider how measurement-based quantum 
computation can proceed using such a state. First, we 
need to verify that single qubit measurements alone suf- 
fice to perform universal quantum computation with an 
encoded cluster state. The logical cluster qubits con- 
sist of four physical qubits, and cluster-state computation 
requires making projective measurements on the logical 
qubits onto entangled superpositions of physical qubits 
the form \0l) ± e~*"|lL). By performing separate pro- 
jective measurements on 3 of the 4 physical qubits in 
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the |±) basis, the logical qubit is decoded into the state 
of the fourth physical qubit with at worst the addition 
of a known Z-eviov if an odd number of |— ) outcomes 
are obtained US ■ Such errors are easily compensated 
for by adapting future measurement bases in the manner 
standard for cluster state computing. 

Next, consider the effect of the perturbative correction 
to the ground state of Eq. ifT^ in terms of its usefulness 
for quantum computation. Clearly, replacing the super- 
position of Eq. H14|l with an incoherent mixture of the 
state |C) and the states \k) can only decrease the useful- 
ness of the ground state for quantum computation. With 
this replacement, the resulting density matrix describing 
the ground state is a mixture of cluster states \C) with 
errors (8) cr^ (and ® cr^) applied to all bonds inde- 
pendently with probability p ~ /{Ag)'^. Note that this 
probability is independent of the size of the lattice, and 
thus these errors can be viewed as independent errors in 
the bonds between sites. 

These errors take the system outside of the logical 
Hilbert space, because a'^ on, say, the first physical qubit 
on a site maps logical states to states spanned by |1000) 
and |0111). However, these errors manifest themselves 
in the single-logical-qubit measurements described above 
as effective Pauli errors in the measurement. Thus, stan- 
dard techniques for fault-tolerance in cluster state quan- 
tum computation |0, 0, Q| can be applied, provided 
the error probability p ~ A^/(4g)^ is below some appro- 
priate threshold. 

Errors can also arise due to the finite energy gap, 
A = 3A^/(16g^). As we are using perturbation theory, 
we require A ^ for the ground state to closely approx- 
imate the encoded cluster state. However, we emphasize 
that the energy scale A need only be small relative to 
g, and not in any absolute sense. The gap A is in turn 
small relative to A. Consideration of the relative magni- 
tudes of these energy scales, in conjunction with deter- 
mining how cold such a lattice can be maintained, will 
determine whether this gap is sufficiently large to allow 
for quantum computation. Specifically, consider the ef- 
fects of using a finite-temperature thermal state of this 
Hamiltonian. As investigated by the thermal state 
will be a mixture of cluster states with logical Z-errors 
occurring independently at each site with probability 
p— (H-exp(A/fc_Br))^^. The critical temperature above 
which the cluster state becomes too noisy to be useful 
for quantum computation scales as fc^Tcrit ~ A |3. 
which in this case is fesTdit ^ {X/g)^X- Thus, the critical 
temperature is determined by the energy scale A and the 
order of perturbation theory at which the degeneracy is 
broken. (For the square lattice, the latter is four, leading 
to the {X/gY dependence.) 

In relation to this point, we note that similar tech- 
niques to those presented here can be used to construct 
systems for any type of cluster state, not just on a square 
lattice, and a large number of other quantum states on 



a graph. This result follows directly from the origins of 
our method in the PEPS formalism. To generalize the 
above method, the number of physical qubits at each 
site must equal the number of bonds, i.e., the number of 
other sites that are directly connected to that site. For 
example, line- or ring-clusters can be created using only 
two physical qubits per site, a hexagonal lattice cluster 
state using three per site, and a cubic lattice cluster state 
using six per site. We note that the number of qubits 
per site determines the order in perturbation theory at 
which the ground-state degeneracy is broken. Continuing 
our argument from the previous paragraph, it is there- 
fore interesting to consider the use of a hexagonal lattice 
cluster state of this form. Such a cluster state is a uni- 
versal resource for quantum computation , and as the 
degeneracy is broken at third order, the energy gap (and 
thus the critical temperature) behaves as Ahox ~ X^/g"^. 
This larger energy gap, as compared with the square lat- 
tice scaling of Agq ~ X'^/g^, may make this state easier 
to prepare and maintain. 

It is worthwhile to consider whether alternative lo- 
cal two-body Hamiltonians also yield an encoded cluster 
state as their ground state, or yield a non-cluster state 
that is nonetheless universal for quantum computation. 
For example, using a Heisenberg antiferromagnetic cou- 
pling on bonds and alternating sites between ZZ-type 
and XX-type Ising interactions leads to a similar result 
as above, where the ground state approximates an en- 
coded cluster state. (To obtain this result, the choice of 
logical basis at each site must be alternated.) In the dy- 
namical approach to cluster state creation, it has been 
shown that the Heisenberg interaction can be used in- 
stead of the Ising interaction ,17.] , and lo gica l encodings 
for this purpose have been investigated [l^l- A future 
line of research would be to investigate if lattices with 
entirely Heisenberg interactions (or other "natural" inter- 
actions on sites and bonds) can yield a computationally- 
useful state as the ground state. We note that our lattice 
of Fig. n and the hexagonal lattice with three physical 
qubits per site, have identical structure to the "CaVO- 
type" lattices and star lattices, respectively, for which the 
Heisenberg antiferromagnetic coupling can lead to exotic 
quantum states [Sl llSl l2i| • 
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